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Abstract 

We present a high precision Monte Carlo study of the finite temperature 
Z<i gauge theory in 2 + 1 dimensions. The duality with the 3D Ising spin 
model allows us to use powerful cluster algorithms for the simulations. For 
temporal extensions of up to N± = 16 we obtain the inverse critical temper- 
ature with a statistical accuracy comparable with the most accurate results 
for the bulk phase transition of the 3D Ising model. We discuss the predic- 
tions of T. W. Capehart and M.E. Fisher for the dimensional crossover from 
2 to 3 dimensions. Our precise data for the critical exponents and critical 
amplitudes confirm the Svetitsky-Yaffe conjecture. We find deviations from 
Olesen's prediction for the critical temperature of about 20 %. 
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1 Introduction 



One of the most interesting predictions of QCD is the existence of a deconfmement 
transition at some critical temperature T c . Finding a precise characterization of this 
transition is, however, still an open problem. In particular it would be interesting to 
have some indication on the order of the transition and, in case it is of second order, 
on its critical indices. The natural framework to pose these questions is that of the 
finite temperature Lattice Gauge Theories (LGT). In this framework, a seminal 
contribution was given more than ten years ago by B. Svetitsky and L.Yaffe 0. 
They showed that, if the deconfining transition of a given (d+1) dimensional gauge 
theory is of second order, then its universality class should coincide with that of the 
ci-dimensional spin model, with the center of the original gauge group as symmetry 
group. This result is usually known as the "Svetitsky- Yaffe (SY) conjecture" and 
has been confirmed in these last ten years by several Montecarlo simulations. Most 
of these tests were made on LGT with continuous gauge groups (in particular SU(2) 
and SU(3)) which are obviously the most interesting from a phenomenological point 
of view, but do not allow in general to reach the high statistics and precision which 
is needed to test the critical behaviour of a statistical system near a second order 
phase transition. To make such a test we follow the opposite route, namely to choose 
the simplest possible nontrivial LGT, and then try to reach in this simplified model 
the highest possible confidence. Therefore we shall concentrate in the present paper 
on the three dimensional gauge Ising model. Besides its simplicity there is another 
important reason to choose exactly the Ising model to test the SY conjecture. The 
only claim of a failure of the SY conjecture that we are aware of was indeed 
obtained in a high precision study of the 3d gauge Ising model with temporal extent 
N t = 8. Such a result, if confirmed, would obviously imply a deep change in our 
understanding of the deconfmement transition. In this paper we shall address this 
question, pushing the analysis up to N t = 16, again in the 3d Ising gauge model, 
using a new original algorithm and very high statistics. We shall show that all our 
results are in perfect and complete agreement with the SY conjecture, we shall also 
give some argument to justify the apparent disagreement found in 0. In order to 
test the SY conjecture we need, first of all, a very precise estimates of the critical 
deconfmement temperature. With this estimate we can test another interesting 
conjecture, due to Olesen 0, which is completely independent from the SY one, 
and is based on the string (or flux tube) description of the infrared regime of QCD. 

A seemingly independent, interesting problem, which is actually related to the 
above one, is that of the so called "dimensional crossover" in layered spin systems. 
Indeed it is well known that the dimensionality of a system as it approaches the 
critical point (and hence the universality class to which the system belongs) is only 
determined by the number of spatial dimensions in which the system has an infinite 
extension. For instance, the Ising model on a lattice with a finite thickness Nt in one 
direction (and the lattice size in the other two directions, say N s , going to infinity) 
belongs to the same universality class of the ordinary two-dimensional Ising model. 
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As Nt increases both the critical indices of the two and of the three-dimensional 
Ising model can be observed. In particular, it is possible to define so called "effective 
critical indices" which smoothly interpolate between the d — 2 to the d = 3 value 
of the indices as the adimensional ratio N t /^(/3) (where £(/?) is the bulk correlation 
length of the system) goes from zero to infinity. The theory of dimensional crossover 
was developed by T. W. Capehart and M. E. Fisher [[| as a particular case of 
the more general theory of finite size scaling theory. Their predictions were then 
compared with high temperature expansions finding in general a good agreement. 
It is interesting to notice that even if ref. H and ref. |lj use different languages and 
are interested in different systems they follow actually the same type of arguments 
(scaling and renormalization group) and reach in fact exactly the same conclusions. 
This is particularly evident in the case of the three-dimensional Ising model, where 
dimensional cross-over and deconfinement transition are in fact the two faces of the 
same phenomenon. Duality, which maps the Ising gauge model onto the spin model 
also relates deconfinement and dimensional cross-over. There is however in H a very 
puzzling result. The critical index A with which the critical temperature on lattices 
of finite thickness approaches the bulk critical temperature turns out to be different 
for different boundary conditions. In particular, in the case of periodic boundary 
conditions in which we are interested here, it does not coincide with the value 
A = \ jv (where v ~ 0.63, is the critical index of the correlation length) which one 
would expect. This is again a result which, if confirmed, would change our approach 
to the analysis of the deconfinement transition, since it is a common and general 
use to extrapolate the finite lattice determinations of the critical temperature to the 
continuum limit, following the scaling behaviour of the correlation length. In this 
paper we shall also address this question, showing that the correct answer is indeed 
A = \ jv and that the result of ||] is probably an artifact of the strong coupling 
expansion. 

A further reason of interest in the 3d Ising model on layered lattices is that in 
these last years this model (and in particular its realization on lattices with very 
few layers) has become a laboratory to test various theoretical methods || || |7j to 
find the critical temperature and the critical indices. Due to the high precision of 
our results we are in the position to test the reliability of these methods and, at 
least in one case (the MHZ approach |7|] , see below) to better understand the reason 
of its failure. 

This paper is organized as follows. Sect. 2 will be devoted to a brief introduction 
to finite temperature LGT and to the Svetitsky-Yaffe conjecture. In sect. 3 we 
discuss dimensional crossover following [|J, while in sect. 4 a brief account of the 
string model approach to the infrared regime of LGT, and the consequent Olesen 
conjecture for the deconfinement temperature will be given. In sect. 5 we shall 
describe the model and the algorithm that we use. Finally section 6 will be devoted 
to a discussion of the results of the simulation and sect. 7 to the comparison with 
other existing estimates of the critical temperature and of the critical indices. 
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2 Finite Temperature LGT 



2.1 General setting 

Let us consider a pure gauge theory with gauge group G, defined on a d + 1 di- 
mensional cubic lattice. In order to describe a finite temperature LGT, we have 
to impose periodic boundary conditions in one direction (which we shall call from 
now on "time-like" direction), while the boundary conditions in the other d direc- 
tions (which we shall call "space-like") can be chosen freely. We take a lattice of 
N t (JV S ) spacings in the time (space) direction, and we work with the pure gauge 
theory, containing only gauge fields described by the link variables <? n;(U G G, where 
n = (x, t) denotes the space-time position of the link and fi its direction. We shall 
denote the time-like direction with the index /i = and the remaining space-like 
directions with roman indices fj, — i — 1 • • • d. We choose for simplicity the same 
bare coupling (3 in the time and space directions. The Wilson action is then 

s w = J2 I P ( E iwGmO + E iwG&w) ] \ , (i) 

n y \ i i<j J ) 

where Tr f denotes the trace in the fundamental representation and g n;0i (g n -,ij) are 
the time-like (space-like) plaquette variables, defined as usual by 

(2) 

In this framework the temperature T is given by: 

where a is the lattice spacing. 

In a finite temperature discretization it is possible to define gauge invariant 
observables which are topologically non-trivial, as a consequence of the periodic 
boundary conditions in the time direction. The simplest choice is the Polyakov 
loop, defined in terms of link variables as 

N t 

P(f) = Tr f y(£)=Tr f n(<fe )t; o) ■ (4) 

t=i 

A second important consequence of the periodic boundary conditions in the 
time direction is that the theory has a new global symmetry (unrelated to the 
gauge symmetry), with the center C of the original gauge group as symmetry group 
(hence, if G is abelian, C will coincide with G). The Polyakov loop is a natural 
order parameter for this symmetry. 

In d > 1, finite temperature gauge theories admit a deconfinement transition 
at T = T c , separating the high temperature, deconfined, phase (T > T c ) from 
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the low temperature, confining domain (T < T c ). The high temperature regime is 
characterized by the breaking of the global symmetry with respect to the center of 
the group. In this phase the Polyakov loop has a non-zero expectation value, and 
it is an element of the center of the gauge group. 



2.2 Svetitsky-Yaffe conjecture 

The idea on which the SY conjecture is based is that, if one would be able to in- 
tegrate out all the gauge degrees of freedom of the original (d + l)-dimensional 
model except those related to the Polyakov loops then the resulting effective the- 
ory for the Polyakov loops would be a d- dimensional spin system with symmetry 
group C. The deconfmement transition of the original model would become the 
order-disorder transition of the effective spin system. This effective theory would 
obviously have very complicate interactions, but Svetitsky and Yaffe were able to 
argue that all these interactions should be short ranged. As a consequence, if the 
transition point of the effective spin system is of second order, near this critical 
point, where the correlation length becomes infinite, the precise form of the short 
ranged interactions should not be important, and the universality class of the de- 
confinement transition should coincide with that of the simple spin model with only 
nearest neighbour interactions and the same global symmetry group. As a conse- 
quence all the critical indices describing the two transitions and all the adimensional 
ratios of scaling quantities should coincide in the limit. This is the case for instance 
of the ratio of partition functions which will be discussed in sect. 5 and of the indices 
v and rj which will be discussed in sect. 6. 

This argument holds also for abelian gauge groups. In particular, the deconfme- 
ment transition of the (2+1) dimensional gauge Ising model which we shall study 
in the following should belong to the same universality class as the two-dimensional 
spin Ising model. 



2.3 Duality 

A peculiar, interesting, feature of abelian LGT in 2+1 dimensions is that they can 
be mapped through duality, into (2+1) dimensional spin models, with the same 
symmetry group (for a review see for instance [§]). 

In the case of the Ising model (in which g n ^ e {±1}) this is accomplished 
by the well known Kramers- Wannier duality transformation which relates the two 
partition functions: 



OC Z spin (P) (5) 

1 

2 



(3 = -J log [tanh(/3)l , (6) 
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where Z gauge (j3) is defined, according to eq.(|I|) and (0) as 

Zgauge(P) = ^ eX P ( X Pfa/w , (<0 

and 

Z apin (P) = E exp(-PH(s)) , (8) 

s*=±l 

with 

H(s) = -J2J(mSiSj (9) 

(«> 

where Sj G {±1}, i and j denote nodes of the dual lattice and the sum is extended to 
the links (ij) connecting the nearest-neighbour sites. For the moment the couplings 
J(ij\ are all chosen equal to +1 . 

Using the duality transformation it is possible to build up a one-to-one mapping 
of physical observables of the gauge system onto the corresponding spin quantities. 
In particular the deconfmement transition in the gauge Ising model becomes the 
order-disorder transition in the spin Ising model. 

It is important to notice that in general the duality transformation does not 
conserve specific choices of the boundary conditions (be in the following). In fact, 
in obtaining eq.(^) (see for instance ||) one usually assumes to be in the thermo- 
dynamic limit and neglects be. In our discussion it is crucial to impose periodic 
be in the time direction in the gauge model, but there is no reason to expect that 
these hold also for the dual spin model. In fact, it can be shown that the corre- 
sponding dual model is then a spin Ising model with fluctuating be in the time 
direction, namely a mixture of the two partition functions with periodic and an- 
tiperiodic be respectively. However the antiperiodic be in the time direction force 
the formation of a space-like interface in the Ising spin model. As a consequence 
the antiperiodic free energy will be depressed with respect to the periodic one by a 
factor proportional to exp(—aN x N y ) and in the thermodynamic limit N x , N y — > oo 
it will become negligible. Hence in this limit the duality transformation actually 
becomes a mapping between a spin and a gauge Ising model, both with periodic be 
in the time direction. This will allow us to make the simulations directly in the spin 
Ising model with periodic be (which is a very convenient choice, since for this model 
(see sect. 5) we have a very powerful algorithm to extract the critical couplings) and 
then to translate the result in the framework of the gauge Ising model, so as to test 
the SY and Olesen conjectures. 

Since this equivalence is of central importance in our work it could be useful 
to discuss it from a different (but equivalent) point of view. We can immediately 
see that a mapping between gauge and spin model both with periodic be is impos- 
sible because there is no dual counterpart for the expectation value of an isolated 
Polyakov loop. Since the Polyakov loop is the order parameter of the deconfme- 
ment transition this seems to exclude the possibility to study this phase transition 
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in the dual model. However the dual of the correlator between two Polyakov loops 
is perfectly well defined. Given any surface in the gauge model bordered by the two 
Polyakov loops the dual of the correlator is given by ZfjZ sp i n where Zf is the par- 
tition function of the spin model in which the links orthogonal to the surface have 
been frustrated. If we now take the thermodynamic limit in the spacelike directions 
(and only in this limit) we can push away as far as we want the two Polyakov loops 
and obtain, by using the standard factorization theorems, the expectation value of 
the isolated Polyakov loop. 



3 Dimensional cross-over in spin models 

In the study of dimensional cross-over in spin models the N t space-like planes are 
usually called layers. Let us assume that the pure d-dimensional model (one layer) 
has a critical point located at f3 = (3 c (N t = 1). If the d spatial dimensions are sent 
to infinity, keeping the number of layers finite, the system still has a critical point, 
whose position (3 = /3 c (N t ) changes continuously with N t . The limit lim^-^oo Pc(^t) 
will coincide with the (bulk) critical point (which we shall simply denote as f3 c in 
the following) of the (d + 1) model. The reduced critical point shift is defined by: 

m) = (10) 

Pc 

We can introduce the critical point shift exponent as 

e(N t ) ~ bi N t - X N t ^ oo (11) 

with b\ a suitable constant. Let us introduce the reduced temperature (with respect 
to the bulk transition of the d + 1 model) : 

Then the bulk correlation length has the following scaling behaviour: 

m ~ h r v N t - oo (13) 

If we assume that the critical point in the layered system is reached when 
£(Pc(N t )) ~ N t then it immediately follows that A = 1/v [|J. Hence we expect 
the value A -1 ~ 0.63 for the cross-over from the three-dimensional to the two- 
dimensional Ising model in which we are interested. 

In view of this analysis the result suggested in [f|] (based on the results of a 
strong coupling analysis) that A -1 ~ 0.63 only for fixed be (in the Nt direction), 
and that A -1 ~ 0.5 for periodic be seems indeed rather surprising. It is in fact very 
unlikely that a change in the be could affect the above analysis, which is essentially 
a bulk result. The only known examples of such a change refer to rather degenerate 
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cases like the spherical model or the crossover from the two-dimensional to the one- 
dimensional Ising model where the choice of the be can indeed drastically affect the 
bulk of the model, changing for instance its operator content. 

Since our measurements have been made exactly in the case of periodic be we 
are in the position to test the result of || . We shall give below strong evidence that 
also for periodic be one has A -1 ~ 0.63. Most probably the result reported in ||J] 
was only an artifact of the strong coupling expansion. 



4 The Olesen conjecture 

It is widely accepted that the infrared regime of a LGT can be described by some 
sort of effective string theory. This conjecture is a natural consequence of the 
observation that the colour flux joining a pair of quarks in the confining phase is 
concentrated inside a tube of finite thickness. According to the conjecture, when the 
quark and the antiquark are pulled very far apart, the flux tube should behave like 
a vibrating string || . Due to the well known problems related to the Weyl anomaly 
even the simplest possible string theory, the bosonic string, cannot be consistently 
quantized in more than one dimension. Depending on the quantization method, one 
finds either the breaking of rotational invariance or the appearance of interacting 



longitudinal modes (Liouville fields). However it is possible to show |I0| that for 
large interquark separation these anomalous contributions can be neglected. The 
resulting effective string theory becomes in this limit an ordinary two dimensional 
conformal field theory. In particular, in this infrared limit, the bosonic string flows 
toward the ordinary (two-dimensional) free gaussian model. In the framework of 
the lattice regularization the interquark potential is given by the expectation value 
of Wilson loops. Thus, according to the above conjecture, we expect the quantum 
fluctuations of large Wilson loops to be well described by some suitable effective 
string model (not necessarily the bosonic one). It can be shown that the contribution 
to the expectation value of the Wilson loop due to these quantum fluctuations can 
be written as a power series in the adimensional parameter (aA)^ 1 where a is the 
ordinary zero temperature string tension and A is the Area of the Wilson loop. 
By using the effective string model one can then evaluate (at least in principle) 
these contributions order by order in the expansion parameter. This approach was 
pioneered in Q I 1 1 (where it was calculated the first term of the expansion) and 
its consequences have been successfully tested in various models by several groups 
during the past years. 

This same reasoning can be extended also to the case of the correlator of two 
Polyakov loops: < P{x)P{x + R) >. From the string point of view the only 
difference with respect to the Wilson loop case is in the be, which can be easily 
taken into account. From the gauge theory point of view we have a new important 
feature: the natural expansion parameter of the theory is now T j ^fo. One can thus 
hope to obtain an universal relation between a and T c by requiring, for instance, 
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that the temperature dependent string tension cr(T) vanishes at the deconfinement 
temperature T c . 



This program was initiated by O.Alvarez and R.Pisarski |12|| , who worked in 
the d — > oo limit (where d are the space-time dimensions). Their result was later 
extended by P.Olesen to any value of d ||. Assuming the bosonic string model as 
effective description, they obtained in the R —>■ oo limit 

< P(x)P(x + R) >~ e - a( - T)R/T (14) 

where cr(T) is given by: 

cr(T) = a* 



T 



(15) 



and 



T > = w^r (16) 

We shall denote in the following eq.([Tj]) as the "Olesen conjecture". 

Let us stress again that this result has been obtained assuming the bosonic string 
as effective model, which could well be a too simple description for the infrared 
regime of realistic lattice gauge theories. For instance in eq. flTUD there is no residual 
information of the gauge group of the model. The only remaining dependence is on 
the number of space-time dimensions of the theory. However, recently, exactly in the 
case of the 3d Ising model, we obtained an impressive confirmation of this bosonic 
string picture fl3|]. In order to describe this result let us first make an important 
observation (see [n| for the details). The bosonic effective string approach described 
above is completely equivalent to the so called "capillary wave model" (CWM) [14], 



which describes the behaviour of fluid (or rough) interfaces in statistical systems. 
This equivalence can be easily understood in the context of the Ising model. The 
same duality transformation which relates the gauge and the spin models, maps 
Wilson loops of the gauge model into interfaces in the spin model (modulo the 
obvious changes in the be), and maps the effective bosonic string into the CWM. In 
fact the CWM, in its simplest formulation, describes the interface with an effective 
Hamiltonian proportional to its area [], and this is exactly one of the possible ways 
to define the bosonic string. In [|13|] we succeeded in evaluating the second term of 
the expansion in a A in the CWM and found a perfect agreement with the results of 
a set of high precision simulations. If this agreement would be conserved to higher 
order of the expansion it would necessarily imply the correctness of the Olesen 
conjecture. This was one of the initial motivations for the present work. Testing 
the Olesen conjecture we are in some sense testing the bosonic string picture to all 
orders (at least in the Ising model). 

1 More precisely, it is proportional to the variation of the surface area with respect to the 
classical solution, which in the Ising case is the cristallographic plane. 
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5 The Numerical Method 



5.1 The model 

We study the spin Ising model denned by eq.@ and (D,ona simple cubic lattice 
with extension N x , N y and N t in the three directions x, y, t. N t will be always smaller 
than N x and N y , and the be in the t and x direction will always be periodic. In the x 
direction we shall implement both periodic and antiperiodic be. Antiperiodic be will 
be obtained by setting the couplings J/^\ = — 1 in the links joining the lowermost 
and uppermost layers in the x direction, while keeping all the other couplings equal 
to +1. In the following we shall denote the Hamiltonian constructed in this way 
as H(s,ap), while the Hamiltonian defined as in eq.(|9|) with all the Juj\ = 1 will 
be denoted by H(s,p). In the following an important role will be played by the 
ratio of antiperiodic and periodic partition functions Z ap /Z p . The reason of this is 
twofold: first, as it was shown in it can be used, as the Binder cumulant, to 
estimate the critical temperature of the system and to extract the critical index v. 
Moreover it was shown in [|L5[| , that this ratio is more efficient and converges faster 
than the standard Binder cumulant. Second, in the present case, owing to the exact 
solution of the two dimensional Ising model, and according to the discussion of the 
previous sections, we can predict exactly the value of the ratio Z ap jZ p at the critical 
temperature. This will give us an independent, very precise method to extract the 
critical temperature. 

In order to obtain the ratio Z ap /Z p one must generalize the system and treat also 
the be as dynamical variables. The partition function of this generalized system is 
given by 

Z = J2J2 exp[-PH(s,bc)] (17) 

be Sj=±l 

where the new index be can take the two possible values: be = {p, ap}. The fraction 
of configurations with antiperiodic be is given by: 



Zap 



{heap) (18) 



and similarly for the periodic ones. Thus we can express the ratio Z ap /Z p as a ratio 
of observables in the same system: 

Z a p (Sbc,ap) / 1n \ 

~v ~ Tk — r> ^ iy J 

£>p \°bc,p/ 

which is hence accessible in a single Monte Carlo simulation, provided the be can 
be updated efficiently This is the case for the algorithm which shall be described 
in the following subsection. 
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5.2 The algorithm 



For the updating of the be we used a specially adapted cluster algorithm which was 
introduced in ref. |16|]. Let us shortly recall the basic steps of this algorithm. 

First the links of the lattice are deleted with the standard probability of the 
Swendsen-Wang [|17] algorithm. 

p d = exp(-/3(l + J^jySiSj)). (20) 

In the second step one checks whether it is possible to change the type of be simul- 
taneously with a sign change of certain spins without changing the sign of J < ij > s i Sj 
for any frozen (i.e. not deleted) link. 

All clusters (in the sense of the standard Swendsen-Wang algorithm) which have 
a site in the x = 1 slice are checked whether they can be uniquely decomposed into 
a part "above" and "below" the x = 1 to x = L boundary. Technically this is 
achieved by introducing an additional label when constructing the cluster using the 
ants in the labyrinth algorithm. The start site of the cluster in the x = 1 slice gets 
the label —1. The sites joining the cluster get the label of their parent site, except 
when the x = 1 to x = L boundary is crossed. Then the new site gets minus the 
label of the parent site. If at the end of the construction of the clusters all sites 
connected by a frozen link have consistent labels the type of the be is altered and 
the spins are multiplied by their label. 



5.3 The Simulations 



We alternated the boundary flip with standard single cluster updates [flql . In the 
case of our largest lattice (128 x 128 x 16) we performed 32 single cluster updates per 
boundary flip. After each boundary flip we performed a measurement of the energy, 
the magnetisation and the type of the be. The total number of measurements was 
100000 throughout. It is interesting to note that the integrated autocorrelation 
time r for the type of the be was smaller then 0.5 due to negative values of the 
autocorrelation function. A similar observation was already made in [plj for the 
L x = L y = L t case. 



6 Discussion 

6.1 Behaviour of Z ap /Z. 



p 

The first test of the Svetitsky-Yaffe conjecture is given by the behaviour of Z ap /Z p . 
In fact, according to the discussion of sect. 2 we expect that at the critical point this 
ratio coincides with the corresponding one in the two-dimensional Ising model. This 
can be evaluated by using the exact solution of the model (see for instance [[0|) or by 
using well known conformal field theory techniques (for a review of this techniques 
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see for instance p0[ ). It depends only on the shape of the two dimensional lattice, 
namely on the ratio between N y and N x . The result is: 

Z^, = ^(0,r) + ^ 3 (0,r)-^ 4 (0,r) 

Z p ^ 2 (0,r)+^ 3 (0,r) + ^4(0,r) 1 ' 

where r = iN y /N x , and ^(O, r) denote the Jacobi theta functions (for the notations 
see for instance ref. ||20|| ). For completeness we remind that Z p is periodic both in 
x and y directions, while Z ap is periodic in the y direction and antiperiodic in the x 
direction. Two interesting limits of this equation are the following: In the N x — > oo 
limit we get ^ 1 as expected, since in this limit the be in the x direction are 

irrelevant. If N x = N y (square lattice) we have = 0.37288488... , a result which 
we shall use in the following. 

In order to test this behaviour we fixed N t = 6, and first simulated the model 
choosing N x = N y = 12,24,48,96. We extracted the critical temperature using 
the usual crossing procedure (see ref. [[L5J) and found the values quoted in tab.l. 
For each crossing we also evaluated the ratio Z ap /Z p (third column of tab.l). All 
the values are in good agreement with the predictions of eq.(pl|). We decided thus 
to assume the S-Y conjecture and extract the critical temperature by requiring 
Z ap /Z p to be equal to 0.37288488.... The results are reported in tab. 2. The gain 
in precision with respect to the usual crossing method is immediately obvious. As 
our best estimate of the critical temperature we take the value corresponding to the 
largest lattice, namely in the present case /3 c (N t = 6) = 0.228818(4). Fixing this 
value of (3 we then studied the system for various asymmetric lattices, and extracted 
the corresponding Z ap /Z p values. The choices of lattice sizes and the results are 
reported in tab. 3 and fig.l together with the predictions obtained from eq. (BTj) . The 
agreement is impressive in the whole range of r values. Up to our knowledge this 
is the first example of such a test of the S-Y conjecture and it is probably the most 
precise and compelling evidence of the validity of the conjecture among those that 
we shall present in this paper. 

6.2 The critical temperature T c 

We then moved to other values of N t keeping the above described approach. We 
always studied lattices with the same values of N x = N y . We determined the critical 
temperature (3 c {N t ) by imposing = 0.37288488.... In tab. 4 we report the set 
of lattices that we studied. For any given N t we chose as our best estimate of the 
critical temperature that corresponding to the largest value of N x . The resulting 
values of j3 c are listed in the third column of tab. 4, while in the fourth column the 
corresponding values for the dual gauge model are given, obtained by using eq.(H). 
We shall postpone to sect. 7 a detailed comparison of our results with other existing 
estimates, but we have already reported in tab. 4 the two sets of estimates: @] and @] 
which we discussed in the introduction. In particular, in the fifth column we have 
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reported the Montecarlo estimates of (for N t = 2, 4, 8) for (3 C . They agree with 
our results, which are however more than one order of magnitude more precise. 
On the contrary the estimates obtained by using high temperature expansions (last 
column of tab. 4) systematically disagree with our results (except for the N t = 2 
case). 

Let us now study the scaling behaviour of these data. Based on general renor- 
malization group arguments we expect corrections to scaling to be of the form 



N t = c x e- v (l + c 2 e A + c 3 e...) (22) 



where A = vuj, with u ~ 0.82 JzlJ] and C\ is related to the constant b\ which appears 
in eq . (|TTD by: c\ = (b\) u . The second correction term reflects the ambiguity in the 
definition of the reduced critical point shift 

m) = TAEAzI, = m(N,)- m = £(JV()(1 _ m (23) 

First we tried to perform a fit with C\, 02,03,1/ and A as free parameters, while 
keeping j3 c = 0.2216546(10) [^T] fixed. However even when taking all results sum- 
marized in table 4 the statistical errors of C2, C3 and A make the result meaningless. 
Therefore we chose a much simpler approach: We solved the equations 

N t = Cl [e(N t )]- v (24) 

and 

N t = d [e'(N t )]- u (25) 

for pairs of N t with respect to C\ and v. The results are summarized in table 5. 
The results obtained from eq.(f2~4"D seem to converge faster than those from eq. (25 ). 



Therefore we regard v = 0.628(3) and l/c% = 1.477(18) obtained from the pair 
Nt = 12 and N t = 16 and using eq. (|24]) as our best approximation. Comparing the 
numbers obtained from eq. (^3) and eq. ( p5|) we still might expect systematic errors 
of the order of 0.005 for v and 0.05 for l/c%. Our best estimate agrees perfectly 
with the existing estimates [y ~ 0.63) of this critical index for the 3d Ising model, 
and disagrees with the estimate {y ~ 0.5) of ref. 0. 

In addition this analysis can be performed using the Zi gauge theory inverse 
temperatures eq.(^) instead of the inverse temperature of the Ising spin model. 

aT c = ^-(p c -Pc{Nt)Y V (26) 
Then we obtain as our best estimates v = 0.6251(26) and aT c = 2.29(3) , to be 



compared, for instance, with the value aT c = 2.2(1) quoted in [22 . 
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6.3 Test of the Olesen Conjecture 



In order to obtain an accurate estimate of the ratio T c j \fa we need the value of the 
string tensionP] at precisely the critical values of 13. A good approximation for these 



values can be obtained taking as reference values those published in [J^J and [|13|] 
and interpolating among them with the law: 

where (3 re f is the coupling at which the reference value of a is taken and f3 is as usual 
the bulk critical point of the 3d Ising model. If \/3 c — /3 re f \ is small enough then the 
uncertainty in the scaling law and in the determination of v can be neglected. We 
can evaluate the error for <j(j3 c ) in a reliable way. The results of this analysis are 
reported in tab. 6. Once the values of <j(/3 c ) are known one can construct the ratio 
T c j y/a and test if it is stable within the scaling region. These ratios are reported 
in the last column of tab. 6. As it can be easily seen scaling is well established. 
Unfortunately this procedure cannot be extended up to N t = 14 and 16, since in 
that region a has never been estimated. As our final result for the ratio we quote 
that corresponding to N t = 12 

T 

1.2216(24) . (28) 



a 



To give an idea of the enhancement in precision with respect to the previous 
estimates we quote the analogous ratio, taken from ref. [p2| : T c j ^fo = 1.17(10). 

Our result disagrees from Olesen's prediction T c /yfa = = 0.977... by 

about 20%. This is not too surprising in view of the remarks made in sect. 4 . On 
the contrary it is rather impressive to see how effective such a simple description 
can be. It would be very interesting to see if some improved string model could 
give better estimates for the critical temperature. 



6.4 Critical indices 

According to the Svetitsky-Yaffe conjecture the critical indices v and i] which char- 
acterize the finite temperature transition point at fixed N t must coincide with 
the corresponding indices of the two dimensional Ising model, whose values are 
v — 1, rj = 0.25. This is a rather important test since in 0, in a set of high 
precision simulations of the 3d Ising model with Nt = 8, a rather larger value for 
7] was found: rj = 0.40(9). This is indeed the only claim for a failure of the S-Y 
conjecture of which we are aware. To extract the two indices we used a standard 

2 As it is well know the string tension of the gauge Ising model and the interface tension of the 
spin Ising model are related by duality and have the same value. To simplify the discussion, in 
the following we shall always refer to a as the string tension. 
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finite size technique, varying the spatial size N x of the lattice. To obtain the index 
v we studied the slope of the Z ap /Z p at the critical point which is expected to scale 
as N^ u . To extract j) we measured the susceptibility x which is expected to scale 
as N%~ v . To extract the indices we compared nearby values of N x . The results 
are reported in tab. 7. We expect that the measured values approach the correct 
asymptotic result as N x increases. It is however interesting to notice that in the 
case of the v index this asymptotic value is approached precociously, indeed all 
the values listed in tab. 7 are compatible within the errors with the expected v — 1 
value. This indicates that, using the Z ap /Z p slope to extract v one has a very small 



correction to the scaling amplitude, a fact that was already noticed in [15]. On the 



contrary in the case of the susceptibility we see a large correction to scaling and 
the expected value is in fact reached only for the largest lattices that we studied. 
Smaller lattices give systematically higher values of t]. This probably explains the 
result of ref. 0. In conclusion we can say that, also for the critical indices, the S-Y 
conjecture is completely confirmed by our high precision test. 



7 Comparison with existing results 

As we mentioned in the introduction, in these last years the layered Ising model 
with low values of N t , in particular with N t = 2, have become a laboratory to test 
theoretical methods @, H, |, 0] and numerical techniques @, |24|, ^5f to estimate the 
critical indices and the critical temperature of statistical systems. In tab. 8 we have 
compared our estimates for the critical temperature with those obtained with the 
mean field transfer matrix approach !J, the Miiller-Hartmann-Zittarz method 0, 
the constrained variational approach and strong coupling expansions 0. Let 
us make few observation on these numbers (for a more detailed description of the 
various approaches we refer to the original literature). 

a] In the transfer matrix version of the mean field approximation || one actually 

deals with clusters of infinite extent. As a consequence the determination of 
the critical temperature is greatly improved with respect to standard mean 
field results^. This explains the impressive precision of the result for N t = 2, 
for which the authors of || were also able to estimate the theoretical uncer- 
tainty. It is very interesting to observe that our result is in perfect agreement 
with this theoretical estimate. As N t increases the method becomes less reli- 
able, however in || the authors give arguments to support the idea that the 
Nt = 3 value is still correct within 4 digits, a result which is again in perfect 
agreement with our estimate. 

b] As it can be seen in tab. 8 the Miiller-Hartmann-Zittarz method (MHZ) which 

was first proposed in ]7| and studied in this context in || gives much worse 
results for the critical temperature. This is a rather puzzling result since 

3 In the case of the 2d Ising model this approach gives indeed the exact result 
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MHZ has some features in common with our approach. MHZ also studies the 
interface free energy which is related to the logarithm of the ratio Z ap /Z p that 
we study. They locate the critical point requiring the vanishing of the interface 
string tension. The main difference between our approach and MHZ is that 
they approximate the Z ap /Z p ratio with the partition function of a suitable 
Solid on Solid model (SOS). In this way they neglect overhangs and handles 
of the interface and exactly this approximation is the reason of the failure as 
N t increases (while it works perfectly in the case of the pure two dimensional 
model, where it gives the exact result). This is indeed very interesting, since 
it gives a further evidence of the relevance of not SOS-like configurations in 



the physics of two dimensional interfaces. This was already noticed in |26 



where the physics of microscopic handles in fluid interfaces was discussed, and 



also in [13] where, in the framework of the Capillary Wave Model (CWM), the 
failure of a simple SOS approximation (which in the CW approach becomes 
a pure free bosonic action) for two dimensional interfaces was shown. 

c] The constrained variational approach || has been applied up to now only to the 
N t = 2 case. The interesting feature of this approach is that it also allows to 
extract the critical indices v and 7. In particular 7 is obtained by looking at 
the variation of the critical temperature as a function of the ratio p between 
the two couplings j3 s and flt- In the region in which we are interested, which 
corresponds to p = 2 in the notations of [EJ the two indices agree with very 
high precision with the values: v = 1, 7 = 7/4 (which imply 77 = 1/4) that 
we have found and that are expected from the SY conjecture. 

As N t increases all the theoretical approaches start to give diverging results. For 
these values of N t the only reliable tool to estimate the critical temperature are the 
Montecarlo simulations. The existing results are reported in tab. 9 and compared 
with our ones. All the estimates quoted in tab. 9 except our ones have been obtained 
by using the method of Binder cumulants. As it can be seen all the estimates agree 
perfectly within the errors. Notice however that our results are in general one order 
of magnitude more precise than the others. This seems to indicate that the method 



that we have used in this work (which was originally proposed in |15|) is indeed 
superior to the standard Binder cumulants method. 
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Tables 



N x 


Pc 


Za.p. / 


12 - 24 


0.22878(5) 


0.376(4) 


24-48 


0.22881(3) 


0.373(4) 


48-96 


0.22882(1) 


0.370(4) 



Tab. 1 Results of the crossing analysis for N t — 6. In the first column the two 
values of N x used to extract j3 c , which is reported in the second column. In the last 
column the value of Z a p /Z p at the crossing point. 



N x 




Slope 


12 


0.22884(4) 


-54.4(4) 


24 


0.22881(2) 


-107.9(7) 


48 


0.228812(9) 


-219.6(1.9) 


96 


0.228818(4) 


-434.1(4.3) 



Tab. 2 (3 C values for N t — 6 and various choices of N x , obtained by imposing 
Za.p./Zp. = 0.37288488.... In the last column is reported the slope of Z a _ p J Z p _ at 

the critical point. 
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N x 


Ny 


Za.p. / Zp_ 


eq-(W 


32 


96 


0.0858(H)(4) 


0.0867 


48 


96 


0.1778(15)(8) 


0.1752 


64 


96 


0.2461(18)(12) 


0.2491 


96 


64 


0.5279(23)(15) 


0.5292 


96 


48 


0.6531(20)(12) 


0.6558 


96 


32 


0.8248(20)(7) 


0.8269 



Tab. 3 Z a . p jZp. values as a function of the lattice shape. In the third column 
the values measured in the MC simulation, the first error is the statistical error, 
the second takes into account the statistical error in the estimate of 
(3 C = 0.228818(4). In the last column the values of eq.(21). 



N t 


N x 




Pc 


i 


1 


2 


4,8,16,32 


0.27604(3) 


0.65608(5) 


0.2758(4) 


0.2760(10) 


3 


24 


0.24607(8) 


0.71102(8) 




0.2449(1) 


4 


8,16,32,64 


0.236025(8) 


0.73107(2) 


0.2358(2) 


0.2347(1) 


5 


40 


0.231421(13) 


0.74057(3) 




0.2301(2) 


6 


12,24,48,96 


0.228818(4) 


0.746035(8) 




0.2273(3) 


7 


52 


0.227195(9) 


0.74947(2) 




0.2258(3) 


8 


16,32,64 


0.226102(5) 


0.75180(1) 


0.2262(2) 


0.2248(1) 


10 


80 


0.224743(5) 


0.75472(1) 






12 


24,48,96 


0.223951(3) 


0.756427(6) 






14 


102 


0.223442(4) 


0.757527(8) 






16 


32,64,128 


0.223101(2) 


0.758266(4) 







Tab. 4 (3 C values for various N t . In the first column the inverse temperature N t , 
in the second column the set of lattices sizes which have been simulated, in the 
third column our best estimates for f3 c (N t ) and in the fourth column the 
corresponding dual values $ c (N t ). In the last two column are reported for 
comparison some existing data for /3 c (N t ) taken from ref j^J (Montecarlo 
simulations) and ref ||/ (strong coupling expansion). 
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N t 


1/ci 


V 


1/ci 


V 


2-4 


1.3043(10) 


0.5903(3) 


1.0393(7) 


0.5208(3) 


4-6 


1.377(3) 


0.6096(7) 


1.230(3) 


0.5824(7) 


5-7 


1.394(11) 


0.6135(23) 


1.276(9) 


0.5935(22) 


6-8 


1.423(8) 


0.6190(16) 


1.323(8) 


0.6035(16) 


7-10 


1.438(14) 


0.6219(24) 


1.357(13) 


0.6103(24) 


8-12 


1.442(10) 


0.6224(16) 


1.375(9) 


0.6134(16) 


10 - 14 


1.438(20) 


0.6218(32) 


1.386(19) 


0.6152(31) 


12 - 16 


1.476(18) 


0.6275(26) 


1.432(18) 


0.6223(26) 



Tab. 5 The critical index v and the constant 1/ci as a function of the pair of N t 
values listed in the first column. In the second and third column the results 



obtained by using eq.(2I). In the last two columns those obtained by using eq.(25) 



N t 


Pc 


2 

a a 




6 


0.228818(4) 


0.0189(2) 


1.212(6) 


7 


0.227195(9) 


0.0138(2) 


1.216(9) 


8 


0.226102(5) 


0.0105(2) 


1.220(12) 


10 


0.224743(5) 


0.00662(12) 


1.228(11) 


12 


0.223951(3) 


0.004654(18) 


1.2216(24) 



Tab. 6 The adimensional ratio T c / \fa for various values of N t . In the second 
column the critical couplings, in the third column the corresponding values of a 
interpolated according to eg. ^\). The reference values of a for the interpolation 
are taken from tab. 2 of ref. for N t = 6, 7, 8 and 10, while for N t = 12 we have 
used the value for (3 = 0.224 reported in j73|/ . 
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N t 


N x 


V 


V 


2 


8-4 


0.979(21) 




2 


16-8 


0.980(21) 




2 


32-16 


1.012(21) 




4 


16-8 


0.971(25) 




4 


32-16 


0.996(21) 




4 


64-32 


1.000(23) 




6 


24-12 


1.012(20) 


0.279(6) 


6 


48 - 24 


0.975(21) 


0.262(6) 


6 


96-48 


1.017(28) 


0.251(5) 


8 


32- 16 


0.979(20) 




8 


64-32 


0.983(21) 




12 


48 - 24 


0.994(24) 




12 


96-48 


0.983(26) 


0.263(7) 


16 


64-32 


0.989(19) 


0.276(8) 


16 


128 - 64 


1.026(26) 


0.263(7) 



Tab. 7 The v and r] indices for various values of N t . In the second column the 
pair of N x values used to extract the index (see text). The susceptibility has been 

evaluated only in a subset of cases. 



N t 


& 


i 





1 


1 


2 


0.27604(3) 


0.276030(2) 


0.2847 


0.2759 


0.2760(10) 


3 


0.24607(8) 


0.2460 


0.2626 




0.2449(1) 


4 


0.236025(8) 


0.2390 


0.2567 




0.2347(1) 


5 


0.231421(13) 


0.2357 






0.2301(2) 


6 


0.228818(4) 








0.2273(3) 


7 


0.227195(9) 








0.2258(3) 


8 


0.226102(5) 








0.2248(1) 



Tab. 8 Comparison of (3 C values obtained with different theoretical approaches 
for small values of N t . In the first column the inverse temperature N t and in the 
second column our best estimates for (3 c (N t ). In the following columns the j3 c 
values obtained with the mean field transfer matrix approach j^j, the 
Miiller-Hartmann- Zittarz method Q], the constrained variational approach j^/ and 

strong coupling expansions j^J. 
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N t 




[24] 


i 


m 


2 


0.27604(3) 




0.2758(4) 


0.27598(1) 


3 


0.24607(8) 


0.24600(8) 






4 


0.236025(8) 


0.23605(5) 


0.2358(2) 




5 


0.231421(13) 


0.23137(5) 






6 


0.228818(4) 


0.22882(3) 






7 


0.227195(9) 


0.22721(4) 






8 


0.226102(5) 




0.2262(2) 





Tab. 9 Comparison of (3 C values obtained with Montecarlo simulations by 
different groups. In the first column the inverse temperature N t and in the second 
column our best estimates for /3 c (N t ). In the following columns the (3 C values 
obtained in ^ and f$BJ respectively. 



Figure Caption 

Fig. 1 Z ap /Z p values as a function of the ratio N y /N x for N t = 6, at the critical 
coupling f3 c = 0.228818. The crosses denote the numerical estimates reported in 
Tab. 3, while the solid line is the theoretical expectation according to eg. (pH[/ 
Errors are not reported since they are smaller than the symbols. 



References 

[1] B. Svetitsky and L.Yaffe,Mfd. Phys.B210 (1982) 423. 
[2] S.Wansleben and J.Zittarz, Nucl. Phys. B280 (1987) 108. 
[3] P.Olesen, Phys. Lett. 160B (1985) 408. 

20 



[4] T.W.Capehart and M.E.Fisher, Phys. Rev. 13B (1976) 5021. 
[5] A.Lipowski and M.Suzuki, Physica A 198 (1993) 227. 

[6] L.Angelini, D.Caroppo, M.Pellicoro and M.Villani, Physica A 219 (1995) 447. 
[7] E.Muller-Hartmann and J.Zittarz, Z. Phys. B 27 (1977) 261. 
[8] R.Savit, Rev. Mod. Phys. 52 (1980) 453. 

[9] H.B.Nielsen and POlesen, Nucl. Phys. B61 (1973) 45; G. 't Hooft, Nucl. Phys. 
B72 (1974) 461. 

[10] POlesen, Phys. Lett. 160B (1985) 144. 

[11] M.Luscher, K.Symanzik and P.Weisz, Nucl. Phys. B173 (1980) 365; 
M.Luscher, Nucl. Phys. B180 (1981) 317 

[12] O.Alvarez and R.D.Pisarski, Phys. Rev. D26 (1982) 3735 

[13] M.Caselle, R.Fiore, F.Gliozzi, M.Hasenbusch, K.Pinn and S.Vinti, Nucl. Phys. 
B432 (1994) 590. 

[14] F.P.Buff, R.A.Lovett and F.H.Stillinger jr., Phys. Rev. Lett. 15 (1965) 621 

[15] M.Hasenbusch, Physica A197 (1993) 423. 

[16] M.Hasenbusch, J.Phys. (Paris) I 3 (1993) 753 

[17] R.H.Swendsen and J.S.Wang, Phys. Rev. Lett. 58 (1987) 86. 

[18] U.Wolff, Phys. Rev. Lett. 62 (1989) 361. 

[19] A.E.Ferdinand and M.E.Fisher, Phys. Rev. 185 (1969) 832. 

[20] C.Itzykson and J-M. Drouffe, Statistical field theory , CUP 1989 

[21] HW.J.Blote, E.Luijten and J.R.Heringa, preprint |cond-mat/95090T6 , 

[22] M.Caselle and F.Gliozzi, Phys. Lett. 273B (1991) 420. 

[23] M.Hasenbusch and K.Pinn, Physica A192 (1993) 342. 

[24] H.Kitatani, T.Kakiuchi, N.Ito and M.Suzuki, Journal of the Phys. Soc. of 
Japan 63 (1994) 2511. 

[25] Z.B.Li, L.Schiilke and B.Zheng, preprint Si-94-09 (University of Siegen) 

[26] M.Caselle, F.Gliozzi, U.Magnea and S.Vinti, Physica A 215 (1995) 21. 



21 



